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PREFACE 


This  report  was  written  to  supplement  material  by 
S.  M,  Pollock  in  Selected  Methods  and  Models  in  Military 
Operations  Research  as  well  as  material  in  OEG  Report  No.  56  by 
B.  O.  Koopman-  Much  af  the  report  is  based  on  a course  given 
by  S . M.  Pollock  at  the  Naval  Postgraduate  School  in  1969. 

The  report's  relation  to  OEG  Report  56  will  be  evident  on 
reference  to  that  work. 

Some  other  sources  of  material  to  which  this  report 
relates  are  listed  in  the  Bibliography. 
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I.  Detection  Theory  and  Detection  Models 

In  signal  detection  theory,  the  decision  making  portion 
of  a detection  system  is  called  the  AcccxueA  and  a detection 
experiment  is  the  observation  by  a receiver  of  its  input  during  a 
time  interval.  The  input  which  is  related  to  a target  is  called 
and  the  input  which  is  not  related  to  the  target  is  called 
naiie..  In  general,  the  observation  is  assumed  to  be  of  a known 
region  which  in  some  cases  is  called  a fizioZutio n celt. 

Signal  detection  theory  is  a basis  for  detection  modeling. 

The  detection  models  which  are  discussed  here  rely  heavily  upon 
it.  When  a detection  experiment  has  been  performed  " ther  the 

event  - {At  least  one  target  was  present  in  the  region  which 

was  observed  during  the  tj.me  of  the  observation.)  or  its  compie- 
ment  Hq  will  have  occurred.  In  general,  the  models  will  specify 
that  either  the  event  Dj^  = (The  receiver  decides  at  least  one 

target  was  present  in  the  region  which  was  observed  during  the 

time  of  the  observation.}  or  its  complement  will  have  occurred. 

In  terms  of  signal  and  noise,  the  events  D^  and  Dj^  can  be 
expressed  as  follows:  = (The  receiver  decides  its  input  during 

the  time  of  the  observation  was  noise.)  and  D^  = {"^he  receiver 
decides  its  input  during  the  time  of  the  observation  was  signal 
and  noise.)  The  theory  which  is  the  basis  for  detection  models 
which  have  the  above  properties  is  called  binary  detection  theory. 

Four  events  which  are  important  in  binary  detection  theory 
are  indicated  in  the  Venn  diagram  of  Figure  1. 
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Figure  1.  Four  events  important  to  binary  detection  theory. 


The  Venn  diagram  emphasizes  a decision  problem  which  is 
associated  with  a receiver.  The  problem  is  this:  Under  what  con- 
ditions should  the  event  occur,  that  is,  under  what  conditions 

should  the  receiver  decide  t.hat  a target  was  present  during  the 
time  of  an  observation?  Detection  theory  may  be  able  to  provide  a 
s-ilution  to  this  problem.  If  this  is  the  case,  the  noise  will 
conerally  be  described  as  a random  process  and  the  signal  will  be 
described  as  either  a deterministic  or  a random  process. 

The  following  notation  and  terminology  will  be  used: 

= P(Dj^JHq)  will  be  called  the  false  alarm  probability. 

• = P(Dj^|H^)  will  be  called  the  detection  probability.  P = P(Hj^) 
.11  be  called  the  prior  probability.  P is  the  probability  that 
target  will  be  in  the  region  to  be  observed  during  the  time  of 
the  observation. 

The  input  to  a receiver  is  assumed  co  be  a quantity  whose 
square  is  proportional  to  power.  The  input  at  some  time  will 

be  sinnbolized  by  y(t^).  The  noise  at  t^  will  be  symbolized  by 
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n(t^)  and  the  signal  by  c(t^).  The  noise  and  signal  processes 
will  be  represented  by  {n(t^),t^et}  and  {s(t^),t^et}  where  t 
represents  the  time  interval  during  which  the  receiver 
observes  the  input.  Often  s(t^)  and  n(t^)  can  be  assumed  to 
be  additive  in  which  case  Hq  = {y(t^)  =n(t^),t^ct}  and 

H,  « {y (t. ) « n(t. ) + s(t. ) ,t.et>. 

11  1 11 

Decause  of  the  finite  quantity  of  information  present  at 
the  input  of  a receiver#  y needs  to  be  measured  at  only  a finite 
number  Ol  points  in  time  in  order  to  be  adequately  determined  over 
an  observation  interval.  For  this  reason#  the  noise  and  signal 
processes  can  be  represented  by  (n(tj^)  » i = and 

{s(t^),  i = l,...#m}  where  t,  are  in  the  observation 

interval  of  length  t.  The  noise  process  is  then  defined  by  a set 
of  m random  variables.  To  specify  the  noise  process,  one  needs 
to  specify  only  the  joint  distribution  of  the  m random  variables 
If  the  signal  process  is  deterministic,  it  is  simply  a finite  set 
of  values  which  is  known  before  the  detection  experiment  is  per- 


formed. 


II. 


Decision  Criteria 


To  simplify  the  discussion  of  decision  criteria  and  decision 
rules,  the  stochastic  process  which  represents  a receiver’s  input 
will  be  assumed  initially  to  be  a single  random  variable  Y.  It 
can  be  called  the  decision  random  variable.  The  input  process  in 
this  case  is  detei-niined  by  the  two  conditional  distribution  func- 
tions FyCylHQ)  and  Fy(y|Hj). 

The  condition  that  a receiver's  input  is  required  to  satisfy 
in  order  that  0^  will  occur  can  be  specified  in  terms  of  a deci- 
sion rule.  For  the  assumed  case,  a decision  rule  is  a rule  which 
determines  for  every  observable  value  of  Y the  decision  that  the 
receiver  is  to  make.  The  decision  rule  can  be  considered  to  be  a 
function  <p(y)  which  relates  the  observable  values  y to  the 
following  two  statements: 

d^:  Decide  the  input  was  noise. 

Decide  the  input  was  signal  and  noise. 

Defining  a decision  rule  <p(y)  is  equivalent  to  defining  a set 
P.  such  that  = (Yen). 

The  problem  which  was  considered  above  can  now  be  restated 
in  the  following  way:  What  criterion  should  be  adopted  in  order 

to  determine  a decision  rule?  A desirable  cliaracteristic  for  a 
criterion  is  suggested  by  the  following  argument:  Consider  the 

ratio  of  posterior  probabilities 

P(H^|y)/P(Hy|y) . 
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One  might  expect  that  the  values  of  y in  should  make  this 

ratio  relatively  large.  But  making  this  ratio  large  is  equiva- 
lent to  making  the  likel'hood  ratio  L(y)  large.  This  suggests 
that  fl  might  be  defined  as  follows:  = (y:  h(y)  ^ k)  where 

K is  some  positive  constant  which  has  yet  to  be  specified. 

Four  specific  decision  criteria  are  defined  next  in  terms 
of  n.  For  each  criterion.  SI  has  the  above  form  with  only  the 
procedure  for  determining  K being  different.  The  decision  cri- 
teria are; 

1.  The  Neyiaan-Pearson  Criterion:  Choose  n so  that  is 

maximized  subject  to  the  constraint  that  s a where  a is 
a specified  value.  For  a continuous  decision  random  variable 


tiic  i.;Oii5tant  K j.9  wnOssii  sO  tiiat  p£  — u. 

2.  The  Bayes  Criterion:  Choose  n so  that  the  expected  cost 

of  a receiver's  decision  is  a minimum.  For  a continuous  decision 


and 


random  variable  K = (C^q-C^q/C^^-C^^)  (1-P) /P  if  > Cqq 

C„,  > C,  , where  C.  • is  the  cost  of  the  event  D.  fl  H • . 

X]  Ij 

3.  The  Ideal  Observer  Criterion:  Choose  n so  that  the 

probability  that  the  receiver  makes  an  incorrect  decision  is  a 
minimum.  For  a continuous  decision  random  variable  K = (1-P)/P. 

4.  The  Minimax  Criterion;  Choose  SI  so  that  the  maximum 
expected  cost  of  a receiver's  decision  is  a minimum.  For  a con- 
tinuous decision  random  variable,  K = (C^q-Cqq/Cq^-C^^) ( 1-P* ) /P* 

if  ^10^  ^00  ^01  ^ ^11'  the  value  of  P which 

would  make  the  expected  cost  of  a receiver’s  decision  a maximum 
if  a Bayes  decision  rule  <p_  were  used. 

B 
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A more  general  discussion  of  the  above  criteria  would  involve 
the  notion  of  a randomizing  rule. 

If  a model  is  adopted  which  specifies  the  conditional  dis- 
tributions and  a decision  rule,  then  the  values  of  and 

c.in  in  principle  be  deterniined.  The  pair  of  values  ^Pf'P(5^ 
called  a receiver  operating  point.  If  the  decision  rule  results 
from  one  of  the  four  likeliliood  ratio  criterion  listed  above,  then 
it  will  involve  the  parameter  K through  the  relation 

fi  = {y:  L(y)  s:K).  And,  for  a given  value  of  K,  since  Q determines 
(p^,p^),  a unique  operating  point  result.  By  varying  K,  a set 
of  operating  points  can  be  generated  whicli  is  called  an  ROC  curve. 
Different  ROC  curves  can  be  produced  by  changing  one  or  both  of 
the  conditional  distributions.  Note,  a change  in  a conditional 
distribution  function  implies  a change  in  the  signal  or  the  noise. 

Decision  ri'l^s  which  result  from  applying  the  four  likelihood 
ratio  criteria  in  a model  in  which  the  input  process  is  deterrained 
by  a set  of  ra  random  variables  can  be  expressed  in  terms  of  a set  ft 
as  follows:  D = ^ ^^1 ' ’ * ' * I-fyj  » * • * ^ where  K is 

specified  in  the  same  way  us  it  is  in  the  corresponding  case  in 
which  the  input  process  is  determined  by  a single  random  variable. 
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III.  Two  Detection  Models 

A particular  detection  model  will  now  be  exeunined.  In 
the  modal,  the  receiver's  input  process  is  defined  by: 


Hq=  {y (t^)  = n (t^) , i = l,...,m}  and 

^1*  ty  - n (t^)  + s (t^)  , i = l,...,m} 

where  the  values  of  y are  measured  every  At  units  of  time. 

The  noise  process  {n(t^),  i = l,...,m}  is  assumed  to  consist 
of  a set  of  m independent  normal  random  variables  each  with 
mean  zero  and  variance  o®.  The  signal  process  {s(t^),  i = l,...,m} 
is  assumed  to  be  deterministic-  Thus,  the  input  process  consists 
of  m independent  normal  random  variables  each  with 

variance  and  each  with  mean  zero  when  the  target  is  not 

present  and  with  mean  = s(t^)  when  the  target  is  present.  Note, 
such  a model  might  be  used  to  obtain  an  optimistic  estimate  of  a 
detection  systems  performance,  since  all  the  information  about 
the  signal  is  assumed  to  be  known. 

The  result  of  the  application  of  a likelihood  c decision 

rule  in  the  model  can  be  expressed  in  terms  of  a rana  ' i.able 

Z.  This  random  variable  is  called  the  cross  correlation  statistic 
and  it  is  defined  by  Z = \ s^Y^.  However,  it  is  more  convenient 
to  express  ;he  results  in  terms  of  a random  variable  V which  is 
defined  by  V = 2/o^.  In  terms  of  this  random  variable. 


and 


p^  = 1 - *(v*) 


Pj  = 1 - ❖(v*-d'i) 
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where  ^ symbolizes  the  standard  normal  distribution  function  and 
V*  = (1/0  ) (o*  In  K + (1/2)  Is.),  dh  ^ I 

The  parameter  d is  called  the  detection  index  and,  since 


Recall,  the  input  stochasric  process  is  assumed  to  repre- 
sent a quantity  whose  square  is  proportional  to  the  input  power. 
Therefore,  the  average  power  input  to  the  receiver  over  the  time 
interval  t during  which  the  receiver  observes  the  input  from  a 
region  can  ba  approximated  by  I y?/m-  The  a^'erage  signal  power 
can  be  approximated  by  S = I st/m  and  the  average  noise  power 
by  I n?/m.  And  the  expected  average  noise  power  can  be  approxi- 
mated by  N = T a^/m  = . In  these  terms,  d = m(S/N). 

The  receiver's  bandwidth  will  be  represented  by  BW. 

Assuming  m corresponds  to  the  value  determined  by  the  sampling 
theorem,  that  is,  m = t/At  = 2t(BW) . the  detection  index  can  be 
written  as  d = 2t(BW)(S/N)  or,  defining  = N/Dw,  as  d = 2t(S/NQ) . 

A model  in  which  the  signal  process  is  not  deterministic 
will  be  considered  next.  The  model  is  defined  by 

Hq:  ^^i’^i'  ^ “ l,...,m> 

^1*  ~ "i  ®i  ' ^ l,...,m}. 

Here  the  noise  process  is  a set  of  m independent  normal  random 
variables  with  mean  zero  and  variance  as  before.  However, 

now  the  signal  process  is  also  a set  of  m independent  normal 
random  variables  with  mean  zero  but  with  variance  The  result 

of  applying  a likelihood  ratio  decision  rule  in  this  model  can  be 
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expressed  in  terms  of  a random  variable  X which  is  defined  by 

X = ^ y?.  When  a target  is  not  present,  the  statistic  X/o*  has 

a chi-square  distribution  with  m degrees  of  freedom  and  when  a 

target  is  present,  the  statistic  X/Ca^+o^)  has  a chi-square 

distribution  with  m degrees  of  freedom.  Hence,  = P(X^^x*/o*) 

and  = P (a^+o* ) } where  is  a chi-square  random 

variable  with  m degrees  of  freedom  and  x*  is  a number  whi ;h 

is  determined  by  the  decision  rule. 

When  the  target  is  not  present,  since  x/a^  has  a chi-square 

distribution,  the  variance  of  X is  2mo'*  and  the  mean  is  mo^. 

When  the  target  is  present,  since  X/(o*+o|)  has  a chi-square 

distribution,  the  variance  of  X is  2m(o^+o^)^  and  the  mean  is 

s 

!n(rr*+cr^),  oy  central  liiuit  theorem,  as  the  number  of  degrees 
of  freedom  m of  a chi-square  random  variable  becomes  large  it 
can  be  approximated  by  a normal  random  variable  with  the  same  mean 
and  the  same  variance-  Hence,  for  sufficiently  large  m, 
and  p^  can  be  approximated  as  follows: 

p^  = 1 - ♦ [ (x*-ma*)/(2ma" ) '] 

and 

“ 1 “ ♦{  t ) - mOg]  / (2mo'* ) 

where  (a^+al)^  has  been  approximated  by  o'*.  This  approximation 
s 

can  also  be  written  as: 

p^  = 1 - ♦ (’/•) 

= 1 - ♦ (v*-d'^) 

with  V*  = (x*-mo^)/(2raa'')  * and  d = (ma^)  ^/2mo'* . 
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Note,  over  the  observation  interval  t,  J s?/m 
the  average  signal  power  so  S = J is 

the  expected  average  signal  power  and,  as  before 
approximately  the  expected  average  noise  power, 
model  the  detection  index  can  be  written  as 

d ~ (S/N)^  * t(BW)(S/N)^. 


is  approximately 
approximately 
N = 0^  is 
Hence,  for  this 
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IV.  Detection  Model  Applications 

When  a likelihood  ratio  decision  rule  was  used  with  either 
of  the  two  ntodels  discussed  above,  the  following  result  was  obtained: 

Pj  * 1 - *(v*) 

= 1 - * (v*-d*^) 

In  the  case  of  a kno%m  signal  (deterministic  process) , a cross 
correlation  receiver  is  required,  and  d = 2t(BW)(S/N).  This  will 
be  called  Case  I.  In  the  case  of  a Gaussian  signal,  a square  law 
receiver  is  required  and,  for  S/N  <<  1 and  (BW)  large, 
d = t(BW)  (S/N)*-  This  will  be  called  Case  II.  In  both  Case  I 
and  Case  II,  the  detection  index  d is  a function  of  (S/N). 

The  two  models  can  be  related  to  radar  and  sonar  systems,  for 
example,  by  using  the  radar  equation  or  sonar  equation  to  deter- 
mine (S/N) . 

In  some  radar  and  sonar  detection  models,  it  is  assumed 
that  the  decision  rule  used  results  in  a required  p^  and  in 
addition  that  a minimum  value  of  p^  is  specified  in  the  sense 
that  a target  is  said  to  be  detectable  only  when  is  greater 

than  or  equal  to  the  minimum  value.  This  minimum  value  of  p^ 
along  with  the  required  p^  define  what  could  be  called  a mini- 
mum acceptable  signal" to-noise  ratio  (S/N)^^^.  In  sonar  models,  this 
minimum  acceptable  signal-to-noise  ratio  determines  the  detection 
threshold  DT  by  the  relation  DT  = 10  log(s/N)^.  If  the  mini- 
mum acceptable  value  of  p^  is  .5,  then  DT  is  usually  called 
the  recognition  differential  RD. 
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A passive  sonar  detection  model  will  now  l>c  discussed.  It 
::>  an  exeunple  of  a model  which  uses  the  concept  of  a minimum  value 
p^.  It  also  illustrates  a method  of  dealing  with  ncn-station- 
. : y noise  and  signal  processes.  In  the  model,  the  event  detection 
the  event  {Xg^s^O}  where  Xg^  represents  signal  excess,  a 
indom  variable  with  expected  value  SE.  The  source  level,  noise 
•vcl,  directivity  index,  and  the  detection  threshold  or  recogni- 
^on  differential  are  also  random  variables  with  expected  values 
NL,  DI,  and  DT  respectively.  The  passive  sonar  equation 
an  be  written  in  terms  of  these  random  variables  or  in  terms  of 
.^ir  expected  values-  In  the  latter  terms,  SE  = FOM  - TL,  where 
1 SL  - (NL-DI)  - DT  is  the  figure  of  merit.  To  determine  the 
: jLaoility  of  a detection  P{Xgg2tOJ,  the  distribution  of  the 
ndora  variable  Xg^.  must  be  specified. 

This  model  could  be  considered  to  be  equivalent  to  one  in 
. ch  the  event  {Xg^tO}  is  the  event  that  is  greater  than 

•'.inirium  accept  able  value.  Altliough  this  definition  may  appear 
e reasonable  than  the  one  which  defines  {Xgj.^0}  as  the  event 
section,  it  would  likely  be  no  more  useful  for  developing  tac- 
s,  designing  exercises,  or  testing  detection  devices. 

The  noise  an-^  signal  processes  should  be  effectively  sta- 
! onary  during  a "look"  at  a resolution  cell,  that  is,  over  an 
Legration  period  if  a Case  1 or  Case  II  model  is  used  to  relate 
, p^,  and  the  signal-to-noise  ratio.  If  either  of  these  models 

used,  the  signal-to-noise  ratio  for  an  integration  period  will 
the  value  of  a random  variable,  and,  therefore,  p^  will  be 
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conditioned  on  the  random  variable  signal -to-noise  ratio.  The 


receiver  is  usually  assumed  to  be  readjusted  so  that 
relatively  constant  from  "look"  to  "look"  at  a given 
cell  in  applications  of  the  model. 


remains 

resolution 
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V.  Search  Theory  and  Search  Models 

Search  theory  deals  with  problems  such  as  that  of  determining 
the  optimal  allocation  oc  search  effort  in  a region  or  that  of 
determining  the  probability  of  detecting  a target  within  a given 
length  of  time  for  a particular  search  plan.  In  treating  these 
problems,  false  alarms  are  often  ignored.  In  effect,  the  time 
spent  on  false  alarms  is  assumed  to  be  iimall  relative  to  the  time 
scale  of  the  problem,  however,  the  cost  associated  with  a false 
alarm  is  not  assumed  to  be  negligible. 

In  general,  it  is  assumed  here  that  during  a search  a 
target  will  move  through  resolution  cells  surveyed  by  the  detec- 
tion system  of  a searcher.  In  the  model  of  such  a search,  the 
probability  that  the  input  from  a resolution  cell  containing  the 
target  will  cause  the  receiver  tc  declare  the  target  is  present 
is  p^.  In  systems  with  opeia. ors,  the  operator  can  delay  the 
decision  regarding  the  presence  or  absence  of  a target  in  a reso- 
lution cell  and  can  recall  the  input  of  adjacent  resolution  cells. 
Such  detection  syster.s  are  difficult  to  model.  To  illustrate 
this,  consider  the  following  hypothetical  system  consisting  of  an 
active  sonar  system  plus  an  operator.  Assume  that  p^  as  a func- 
tion of  the  signal-to-noise  ratio  has  been  determined  for  the 
operator  in  a laboratory  experiment  by  forcing  the  operator  to 
respond  after  a single  "look"  at  a resolution  cell.  Operationally, 
the  probability  that  the  operator  detect,  a target  will  not  be 
given  by  the  laboratory  determined  value  of  p^  if  the  operator 
delays  the  detection  decision  for  several  "locks."  Various  models 
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havs  been  proposed  to  deal  with  this  situation.  In  one,  a "three- 
out-five"  detection  criterion  has  been  adopted.  The  criterion 
states  that  an  operator  will  declare  a target  if,  out  of  five 
consecutive  inputs  from  adjacent  resolution  cells,  at  least  three 
are  such  that  they  would  have  caused  the  operator  to  declare  a 
target  present  in  a single  "look"  experiment.  An  input  which 
satisfies  this  condition  is  referred  to  as  a success.  With  this 
criterion,  the  probability  that  a target  will  be  detected  with 
m consectuive  "looks"  at  the  target  is  found  as  follows:  Deter- 

mine the  2®  sequences  of  successes  and  failures  of  length  m 
which  could  result  with  m consecutive  "looks."  The  probability 
is  equal  to  the  sum  of  the  probilities  of  occurrence  of  those 
sequences  for  which  the  three-out-of-f ive  criterion  is  satisfied. 
The  probability  of  a particular  sequence  will  depend  on  the  value 
of  and  1 ~ P<j  each  of  the  «i  "looks"  in  the  sequence. 

Unfortunately,  in  most  situations,  the  number  of  sequences  to  be 
considered  is  too  large  to  be  computationally  tractable  even  if 
the  probabilities  in  the  sequence  can  be  determined. 


15 


VI.  The  Probability  of  Detection  During  a Search 

The  problem  of  determining  the  probability  that  a target 
will  be  detected  in  a search  at  or  before  n "looks"  have  occurred 
is  generally  basic  to  the  solution  of  search  problems.  Let  N 
represent  the  number  of  the  "look"  at  which  detection  first  occurs, 
then  this  probability  can  be  expressed  as 

n 

P(N*in)  = I P(N=i). 
i=l 


It  can  also  be  expressed  as 


n 

P(N«n)  = 1 - “IT  (l-g  ) 
i=l  ^ 

where  g^  = P(N=i |N*i-l)  is  the  probability  that  the  event  detec- 
tion occurs  on  the  i—  look  conditioned  on  the  event  detection 
has  not  occurred  earlier.  The  second  expression  for  P(Na»n)  is 
generally  of  greater  interest  than  the  first  one,  since  g^  can 
be  directly  related  to  operational  par^unetcrs  such  as  target  range. 
Note,  if  the  resolution  cell  being  examined  on  the  i—  look  is 
empty,  then  g^^  will  be  zero. 

A continuous  analog  to  P(Ni^n)  can  be  developed  as  follows: 
Consider  a model  in  which  the  time  for  a look  is  At,  and  if 
detection  occurs  on  the  i—  "look,"  the  time  of  detection  is 

•J*  Vi 

lAt-  The  probability  of  detection  at  or  before  the  n — "look"  then 
can  be  expressed  in  i.erms  of  time  by  noting  that  P(NSn)  = P(T*nAt) 
where  T represents  the  time  of  detection.  Now  define  the  detec- 
tion rate  to  be  y (nAt)  = (1/At)  P [T=nAt  lT<:(n-l  J At]  . The  probability 
P(T^nAt)  can  be  approximated  in  terms  of  a related  quantity  y(t)  by 
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P(Tiit)  » l-exp{-  I Y(T)dT} 

0 

where  Y(t)  is  related  to  y (n&t)  by  the  requirement  that 
P(T*t)  equals  P(KiCn)  for  t » nAt.  in  some  models  where  it 
would  be  useful  to  use  P(T*t),  the  continuous  approximation  to 
P(N*n)*  it  miqht  not  be  reasonable  to  model  the  detection  process 
in  terms  of  single  "look"  probabilities.  If  this  is  the  case, 

Y (t)  must  be  defined  by  some  other  means.  The  visual  detection 
model  developed  in  the  OEG  Report  No.  56  is  an  example  of  this. 

Generally,  in  search  operations,  the  target  will  be  moving 
on  a relative  track  C.  In  some  cases,  it  may  be  useful  to  con- 
sider the  track  to  consist  of  segments  C,  ,C-.,....C  . Let  F. 
represent  the  event  the  target  is  first  detected  on  track  segment 
C^.  The  probability  of  first  detecting  a target  while  it  moves 
along  C can  then  be  expressed  as  follows; 

n 

P(F)  =:  l-exp{-  I F(C.)> 
i=l  ^ 

where  F(C^)  is  called  the  sighting  potential  on  the  track  seg- 
ment C.  and  exp{-F(C^)}  = P(r^|f^  ...  . It  is  conven- 

ient to  use  the  concept  of  sighting  potential  in  the  analysis  of 
some  search  types,  for  example,  "ladder  searches". 

The  track  of  a target  which  moves  in  a plane  can  be  described 
by  equations  x » x(t)  and  y * y(t).  These  equations  describe 
the  relative  track  of  the  target  with  respect  to  the  searcher 
(often  called  the  observer)  if  x and  y refer  to  a coordinate 
system  in  which  the  searcher  is  located  at  the  origin,  in  this 
case,  r the  range  of  the  target  is  given  by  r=  Ix^+y®]**. 
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. If  in  a model  the  detection  capability  of  a searcher 

against  a target  is  assumed  to  be  independent  of  the  time  at  which 

the  target  occupies  a particular  position  on  the  target's  track, 

then  time  can  be  replaced  by  target  distance  along  the  track. 

ft 

This  can  be  done  by  using  the  transformation  equation  s = w(T)dT 

■>0 

where  w is  the  target's  speed  relative  to  the  searcher  and  the 
zero  of  time  has  been  chosen  to  coincide  with  the  start  of  the 
search.  If,  in  addition,  the  detection  capability  of  the  searcher 
against  the  target  is  assumed  to  depend  only  on  the  range  of  the 
target,  then  the  sighting  potential  for  a target's  track  C for 
a continuous  model  can  be  expressed  in  terms  of  either  t or  s 

rt  rS 

as  follows:  F(C)  = 1 Y(r(T)]dT=  1 Y (r  (c)  J dcr/w  (c)  sxiice 

•'O  •’0 

dt/ds  = l/w(s) . 
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VII.  Experimental  Validation  of  Detection  Models 

Consider  a detection  systeia  which  could  be  satisfactorily 
modeled  by  a continuous  detection  rate  function  that  depends  only 
on  r,  the  target's  range,  if  a suitable  detection  rate  function 
could  be  found.  Suppose  on  some  basis  that  a particular  detec- 
tion rate  function  is  propos<jd?  An  estimate  of  the  suitability 
of  the  proposed  function  could  be  made  using  standard  statistical 
tests  if  appropriate  operational  data  could  be  obtained.  For 
example,  the  number  of  detections  at  various  ranges  for  straight 
line  encounters.  If  such  data  were  available,  its  agreement  with 
the  cumulative  distribution  function  of  the  target's  range  at 
detection  deduced  from  the  proposed  detection  rate  function  could 
be  determined  by  use  of  the  Kolxaogorov-Smirnov  test.  This  technique 
could  be  applied  even  though  the  number  of  encounters  in  which  a 
target  is  not  detected  was  unknown.  In  ttiis  case,  the  cumulative 
distribution  function  to  be  compared  is  Fj^(rjD)  where  D = 

{detection}  and  R represents  the  target's  range  at  detection. 

To  derive  Fp(riD5  or  equivalently  the  density  function 
fj^(r|D)  given  a continuous  detection  rate  function,  one  can  pro- 
ceed as  follows:  First  find  the  joint  density  function  which  de- 

scribes the  distribution  of  the  target's  rectangular  coordinates 
(X,Y)  at  its  detection  point.  To  do  this,  se>,  f^  y(x,y)  lD)AxAy  « 
fy  y(x,y)AxAy/p  where  p = P(D) . This  can  be  done,  since 
[ (Y=y)  fid]  •=  (Y=y)  . That  is,  Y can  be  considered  to  take  on  a 
value  only  when  a detection  occurs.  Note,  p is  the  value  of  the 
integral  of  ^ix,y)  over  all  pairs  of  values  (x,y) . Now  note  that 

A#  1 


“ fy|j^(y|x)  f jj(x)  and  that  the  conditional  density 
fyly{y|x,  can  be  written  as 

fy 

fyix<ylx)  = (l/w) Y[r(y) lexp{-j  Y Ir (n) Idn/w) . 

—CO 

Given  Y(t)«  fj^(r|D)  can  be  found  by  first  trans- 

forming to  r and  o by  using  the  transformation  equations  x = r 
cos  o and  y = r sin  o and  then  finuxng  the  marginal  distribution 
of  R as  follows: 

2it 

0 


f 


X 


An  example  of  the  above  procedure  will  now  be  given. 

/ i/2a  jxj  « a . kli/r*  jxj  » 

(x)  = I and  y{x)  = | 

' 0 |x|  > a I 0 |x|  > 


Suppose 

a 

a 


The  transformation  to  r and  o gives 

r sin  o 

I (l/2awp) (kh/r*)exp{-| 


kh  (r^cos^+n^ ) ^'^^dn/w} 


fR^j:(r,a|D)  = 


r cos  o * a 


r cos  o > a 


For  large  a,  this  can  be  simplified  by  letting  a * after 

f ^ 

noting  that  p = l/2a  p(x)dx  - W/2a  where  W is  the  sweep 
width.  In  tlie  limit,  ^ fj^  j.(r,crlD)  * (kh/wWr“ ) exp{-kh/wr^i-sin  o)}, 
and  after  performing  the  integration  over  o, 

fj^(r[D)  = (W/2r^>  U-»(kh/wr^)*"}. 

Note,  which  emphasizes  the  fact  that  in 

general  detection  need  not  take  place  during  a straight  line  encounter. 
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VIII.  Sweep  Width  Determination  for  a Random  Track  Angle 

Now  consider  the  problem  of  determining  probability  p of 
detecti.ig  a target  during  an  encounter  which  could  be  described 
by  a straight  line  encounter  model  if  the  angle  <p  between  the 
target's  track  and  searcher's  track  were  known.  Suppose  it  can 
be  assumed  that  lateral  range  and  the  relative  speed 
w - [u^+v^-2uv  cos  g]  are  values  of  random  variables  whose  joint 
distribution  is  known.  The  probability  p of  detecting  a target 
during  a straight  line  encounter  can  then  be  expressed  as 

aOi>  .•  OO 

P = j p(x,w)fy  ^,{x,w)dxdw 

_ot.  0 

where  p(x,w)  = p (det ]x  - x n f = w)  and  ^(x,w)  is  the  joint 

distribution  of  the  lateral  range  X and  the  relative  speed  V. 

If  X and  y are  independent  random  variables  and  if 
p(x,w)  as  a function  cf  w is  adequately  described  by  a linear 
approximation,  then  p can  be  approximated  as  follows: 

p = I p (x,w) f^(x)dx 

— OO 

v.’bere 

w = E(^i  = I wf^(w)dw  - (u^+v^-2uv  cos  g)  ^ f^(g)dg. 

•'0  •'0 

An  average  sweep  width  can  be  defined  by 

r 

w - p(x,w)ax. 


— 00 

If  i.s  uniformly  distributed  so  ~ value 

of  an  elliptic  integral  of  tne  second  kind.  An  average  sweep  width 
based  on  this  value  of  w should  be  used  in  the  Random  Search 
Formula  when  either  u or  v are  not  appropriate. 
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IX.  A Parallel  Sweep  Search  Model 

A model  of  i.  parallel  sweep  or  **iadder*'  search  will  be 
considered  next-  In  the  model,  the  target's  velocity  is  zero 
and  its  position  is  uniformly  distributed  in  a rectangular  region. 
The  searcher  moves  along  a series  of  n parallel  tracks  to  that 
n straight  line  encounters  will  occur  in  a complete  search  of 
the  region-  The  track  spacing  is  s and  the  tracks  are  numbered 
from  1 through  n.  The  search  geometry  is  shown  in  Figure  2. 


n 


r~ 

ns 


-♦  z 


Figure  2.  Search  geometry  for  a parallel  sweep  model. 

If  the  target's  horizontal  coordinate  is  a random  variable 

Z in  the  coordinate  system  shown  in  the  figure,  then  P(det|Z=  z) 
n 

= l-exp(-  F(x.)J  where  x.  = is  the  target's 

i = i 3 J 

J th 

lateral  range  from  the  j — track  and 


ns 


P {det I Z=z) dz/ns . 


P (det) 


In  some  cases,  it  is  convenient  to  use  the  following  coordinate 
system:  Let  I be  the  random  vi^iriafale  wnicfc  is  the  number  of 

the  track  closest  to  the  target.  Relabel  the  tracks  so  that  the 
number  of  the  track  cic-se.it  to  the  target  is  P.  then  = -js  + x 
where  now  j - ..,M  where  n = M + K + 1 and  the  value 

of  K and  M are  determined  by  the  value  of  I,  that  is,  by 
the  location  of  the  target.  Now 


P (det 1 X = X 


Atl 

ni=i)  = l-exp[“  ]]  F(x-)]  = p(x  i) 

j=-k  ^ 


Hence, 


rs/2 

P(det|l  = i)  = I p(x,i)dx/s  ~ p(i). 
-s/2 


n n 

P(deM  = I p(i)P(I=i)  = (I/i.)  I p(i) 
i»l  i=l 


n fs/?. 

P(det)  - (1/ns)  ^ j p(x,i.dx 

-s/  2 


If  n is  large  and  the  detoction  law  is  such  that  y(r)  is 
essentially  zero  for  targets  at  rai.ges  equal  to  or  greater  th-n 
one  or  two  track  spacings  away,  then  for  alTiCSt  all  values  of  i 


and  j,  p(i)  wii:*  be  approximately  equal  to  p(j)  and  it  can 

fs/2 

be  assu.aed  that  p(det)  = (N/Ns)  p(x;s)dx  where 

-s/2 


p(x;s)  = l-erpi-  I F(-js+x,]. 

j =1  — QO 


•>  ^ 


To  illustrate  an  application  of  the  above  result f suppose 
Y(r)  = kh/r^.  Then  F(x^)  = 2kh/wx ^ . If  it  is  assumed  that  the 
number  of  tracks  and  the  track  spacing  is  such  that  p(i)  = p(j) 
for  almost  all  values  of  i and  j,  uhen 
m " 

I F(x.)  = F(-js+x)  « (2kh/w)  (ir*/s^)  csc^  <rx/s) 
j=-k  ^ j=0 

so 

p(x;s)  =■  1 - exp{- (2khir^/ws^ ) csc^  (xx/s)  } . 

and 

P(det)  = 24'{(x/2)  ’'^  (W/s)  } - 1 
1. 

where  W = 2(27rkh/v)  * is  the  searcher*  s sweep  width  against  the 
target.  The  ratio  W/s  is  called  the  coverage  factor.  A plot  of 
P(det)  is  given  as  a function  of  W/s  in  the  National  Search 
and  Rescue  Manual. 

The  track  length  which  is  expended  in  an  unsuccessful  search 
is  nb  by  the  assumptions  of  the  parallel  sweep  model.  Since 
the  area  searched  is  nsb,  the  probability  of  detecting  a target 
with  a random  search  in  the  area  with  track  length  nb  is  given 
by  P(det)  = 1 - exp{-W/s}.  This  value  could  be  compared  to  that 
given  by  the  parallel  sweep  model  in  order  to  obtain  a relative 
measure  of  effectiveness  for  a parallel  sweep  search. 
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X . The  Optimal  Allocation  of  Search  Effort  Problem 

The  problem  of  determining  an  optimal  allocation  of  effort 
for  a search  is  discussed  next.  Here/  an  allocation  of  effort  is 
an  optimal  allocation  if  it  maximizes  the  probability  of  detecting 
the  target. 

The  class  of  searches  which  will  be  considered  aie  those 
which  can  be  described  by  a random  search  model  that  satisfies  the 
following  conditions:  The  target  is  fixed  in  a region  of  area  A. 

The  region  consists  of  m subregions  which  are  determined  by  the 
conditions  that  in  the  i~  subregion  the  sweep  width  is  a con- 
stant w^  and  the  probability  density  function  is  a constant 
” ^i'^^^i  where  is  the  area  of  the  subregion  and  is 

the  target's  prior  probability  of  being  there.  (The  subregions 
are  assumed  to  be  numbered  so  that  W,  P,  >W-P,  > . ..  > W P .)  If  a 
target  is  detected  in  a subregion,  the  searcher  must  be  in  the 
subregion  when  it  is  detected.  The  probability  of  detecting  the 
target  is  then 

n 

P(det)  = I (1-exp  ( } Pi 
i=l 

where  is  called  the  search  effort  density. 

Kow  consider  the  problem  of  determining  the  solution  set 
which  maximizes  P(det)  subject  to  the  constraints 
i = \ *nd  <p^  S 0.  Note,  £ is  the  available  track 

length.  This  is  a nonlinear  programming  problem  for  which  the 
solution  set  is: 
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i in  Q 


InW.  p.  -d/X)  I (AVW.)ln  p.  + «/X 

0 i not  in  SI 

where  X = / A./W.  and  SI  = where  k is  determined 

by  the  condition  that  if  k + i were  included  in  Si  then  0* 

Suppose  the  model  described  above  were  used  to  determine 

an  optimal  allocation  of  track  length  for  a search.  The  result 

would  be  a set  of  values  9.^  = (A./W.)4>-  which  maximize 

1 op  i 1 op 

the  probability  of  detecting  the  target.  Given  detection  will 
occur,  in  what  order  should  the  subreg3.ons  be  searched  in  order 
to  minimize  the  expected  time  until  the  target  is  detected?  The 
following  procedure  will  effectively  minimize  the  expected  track 
length  to  detection  and  gives  a suitable  order  assuming  the 
searcher's  speed  remains  constant  so  that  minimizing  the  expected 
track  length  to  detection  is  equivalent  to  minimizing  the  expected 
time  to  detection. 

Divide  S,  into  minimum  increments  A£  consistent  with 
the  random  search  model  and  then  allocate  the  increments  in  the 
following  order:  Allocate  AIL  to  the  Region  1.  This  will  be 

assumed  to  be  effectively  consistent  with  the  prescription  deter- 
mined by  the  formula  for  tp . given  ASL.  If  this  search  is 

1 op 

unsuccessful,  determine  the  optimal  allocation  prescribed  for 
2A8,.  This  will  effectively  indicate  whether  the  next  increment 
AH  should  be  assigned  to  Region  1 or  Region  2.  Continue  in  this 
fashion  until  detection  or  until  H has  been  expended.  Note, 
this  is  an  approximate  procedure  in  the  sense  of  the  model. 
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For  a discussion  of  the  above  procedure,  let  E(L|Laii) 
symbolize  the  expected  track  length  to  detection  given  detection 
will  occur  by  the  time  track  length  i has  been  used.  That  the 
above  search  procedure  will  effectively  minimize  E(LjL*il)  for 
an  optimal  allocation  of  track  length  is  suggested  by  the  follow- 
ing argument;  Divide  1 into  segments  A?,  as  indicated  above, 

then  the  probability  that  the  target  is  detected  on  or  before  the 
• ttl 

1 — step  of  the  search  P(La;iA)c)  is,  in  effect,  a maximum, 

since  no  other  allocation  of  search  effort,  i will  give  a 

significantly  greater  value  for  any  value  of  i.  Hence, 

Fj^(iAi|i>£)  = P (L^^iAft  |l*S.)  is  also  effectively  a maximum  for  any 

value  of  i since  P (L^iAJ,  = P (lAiA  ? ) /P  (L^  • ) and  P(L*i) 

must  be  equal  to  its  maximum  value  for  any  alternate  procedure. 

n 

Therefore.  £(L|LsiA)  = I (1-F.  (iA£  |Lsi)  ) is  effectively  a 

i=l  ^ 

minimum  for  the  procedure. 

The  above  procedure  is  the  one  which  would  have  been  followed 
if  after  each  unsuccessful  search  one  calculated  new  prior  proba- 
bilities using  Bayes  procedure  and  then  did  the  search  prescribed 
for  a track  length  AH  with  the  new  priors.  for  a discussion  of 
this  point,  see  OEG  Report  No.  56.  An  application,  in  effect,  is 
discussed  in  Reference  4, 
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XI.  Target  Position  Probability  Distributions  Which  Change  in  Time 
C....isider  the  following  model  in  which  a target  moves  with 
constant  course  and  speed  in  a plane  with  its  position  specified 
in  a rectangular  coordinate  system:  The  target's  course  and  speed 

are  independent  random  variables  4 and  U respectively 
whose  density  functions  ^nd  are  given.  The  joint 

density  f;  »•  >.on  fj^  y(X/yrt)  of  the  target's  coordinates  is 
needed/  but  ly  f^^  y(x,y;0)  the  joint  density  function  for  time 

zero  is  given. 

For  such  ;ases/  f^  found  in  principle  as 

follows:  Fir*-  . it  is  convenient  to  use  polar  coordinates  and  with 

X = p cos  a and  y - p sin  a and  to  set 

f(p,o;t)  = f y (p  cos  a,p  sin  a;t) . 

Note/  to  first  order,  P{ target  is  in  AA  at  t)  “ f(p,o;t)AA. 


Figure  3.  Problem  geometry. 
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Referring  to  Figure  3, 


f (p,a;t) hA 


J ] f [r{u,*)  ,e(u,4.)  ;0]AA’f^(^)fu(u)d«du 
0 0 


where  the  functions  rCu^^*)  &nd  6(u,0)  are  such  that  if  the 
target's  polar  coordinates  were  r(u,^)  and  6(u,^)  at  time  0, 
then  with  the  speed  and  course  u and  <j>  they  would  be 
(p,a)  at  t.  The  functions  can  be  expressed  in^licitly  by 


r"  = p®  + (ut)*  - 2put  cos(a-(|i)  and  r cos  0 = p cos  a - ut  cos 


The  integration  is,  in  effect,  a sum  over  all  possible  pairs  of 
values  of  u and  Each  pair  determines  a vector  ut  which 

translates  target  positions  at  time  0 to  target  positions  at 
time  t such  that  AA  is  AA’  translated  without  rotation  or 
distortion.  So,  AA*  - AA,  independent  of  u and  Therefore, 


n2n 

f Ir  (u,<j))  ,0  (u,(l))  ; 0]  ((> ) (u ) d(tdu 


0 0 

As  an  example,  suppose 


fx,Y(x,y;0)  = e 


with  f 


♦ 


'/  = 


2¥ 


, 0 < ♦ s 2ir,  but  with  U =*  u known.  In  this 


case. 


r^  (Uj<|>) 

f Ir(u,(J>)  ,0  (u,<fr)  ;0J  * 2^-r  e 


and,  using  r*  = p^  + (ut) ^ - 2p  ut  cos (a-^)  , 
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f (p,a;t) 


J 


2it  - 5^  [p*4  (ut)  -2p  ut  cos  (0-4))  1 


2iro 


T e 


P ut 


1 “ 2^  to  +(»it)  1 1 f2u  ^^cos(oi-*) 


1 ' 5o^  . 

■f-  e Jo  ^“5^^ 


2ttg 


where  is  the  zero  order  Bessel  function  of  the  first  kind 

and  i = To  obtain  fj^  y(x,y;t),  substitute  P*  = x*  + y* 

in  the  above  result.  In  this  case,  since  fj^^Y  symmet- 

ric about  the  origin,  the  marginal  distribution  of  the  target's 
range  R from  the  origin  is  of  interest,  in  general. 


fR  0(P,a;fc)  = Pfjj  Y^P  ® Pf(P.a;t) 


and 


f2^ 

J 
0 


fR(P;t)  = J €r  g(P,a;t)da 


so,  in  this  case. 


= 5-*-  ® 


P - 2F- 


The  distribution  function  FrCp;!)  for  this  case  is  given  in 
Reference  5. 
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